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@ Signaltrennungsverfahren und -einrichtung fur nichtlineare Mischungen unbekannter Signale 

(57) Die Erfindung ist ein Verfahren und eine Einrichtung zur 
Separierung eines Signalgemisches, bestehend aus der 
nichtfinearen Oberiagerung von M unbekannten Signalen, 
dadurch gekennzeichnet, das N (N ^ M) unterschiedliche 
Mischungen der M Quellsignale inklusive eines eventuelien 
Storsignals einer Signalauswerteeinrichtung zugefuhrt wer- 
den, die durch eine statistische Analyse der Signale die 
nichtiinearen Ubertragungsfaktoren bestimmt und mix die- 
sen errechneten Faktoren die Mischung ruckgangig macht, 
so dalS die N Ausgange der Signaitrennungseinrichtung die 
M Quellsignale ohne Oberlagerungen enthalten. Dadurch 
wird erstrnals die Behandlung nichtlinearer Gemische mog- 
lich, wobei nichtlinear bedeutet, da& die Quellsignale durch 
etn unbekanntes nichtlineares System gemischt werden. Das 
unbekannte System wird durch eine Volterra-Reihe be- 
■ schrieben, und die Signaitrennungseinrichtung bestimmt die 
k Koeffizienten der Volterra-Reihe. Mit derer Kenntnis ist eine 
Entrnischung moglich, au&erdem konnen die Koeffizienten in 
> weiteren Analysen zur Orts- oder Geschwindigkeitsbestim- 
mung der Signalquellen benutzt werden. 
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Beschreibung 



1. Technisches Gebiet 

In der MeBtechnik besteht fast immer das Problem, daB bei einer Messung nicht nur Signale einer einzelnen 
Quelle, sondern oft nur Mischungen aus den Signalen mehrerer Quellen durch die Detektoren aufgenommen 
werden konnen. 

jo 2. Stand der Technik 

Eine Signaltrennung mit moglichst wenig a-priori-Information liber die Signale wurde von Herault und Jutten 
in den folgenden Veroffentlichungen behandelt: C Jutten and J. Herault, "Blind separation of sources, part I: An 
adaptive algorithm based on neuromatic architecture", Signal Processing, vol 24, No. 1 pp 1 — 10, July 1991; 

is P. Common, Q Jutten, J. Herault, "Blind Separation of sources, part II: Problems statement" Signal Processing, 
vol 24, No. 1, pp 11—20, July 1991, E. Sorouchyari, C. Jutten, J. Herault, "Blind Separation of Sources, part III: 
Stability analysis, ", Signal Processing, vol 24, No. 1, pp 21—29, July 1991. Dabei wird das Problem von mit 
konstanten Faktoren linear uberlagerten, zeitlich unverzdgerten, unabhangigen Signalen beschrieben. Wie sich 
herausstellt, funktioniert das Verfahren nicht mit GauBverteilten Signalen, das sind zum Beispiel Sprache oder 

20 Musik. 

In der Veroffentlichung von Molgedey und Schuster (Molgedey L. und Schuster H. G. 1994: Separation of a 
Mixture of Independent Signals using Time Delayed Corelations, Phys. Rev. Lett. 72, 3634—3638.) wird die 
Trennung fur unkorrelierte Signale, die zeitlich unverzogert linear uberlagert sind, behandelt, die auch fur 
GauB'sche Quellen benutzt werden kann. AuBerdem wird ein Verfahren zur Trennung von Mischung skizziert, 

25 bei denen die Quellsignale (zeitlich unverzogert) multipliziert sind. 

Gilles Burell (Burell, Gilles 1992: Blind Separation of Sources: A Nonlinear Neural Algorithm, Neural Net- 
works, Vol. 5, pp 937—.) schlagt einen weiteren Ansatz fur die Trennung von Gemischen aus unabhangigen 
Signalen vor, in dem ein neuronales Netz speziell fiir eine bekannte Art der nichtlinearen, zeitlich unverzdgerten 
Mischungen konstruiert wird Das Verfahren zum Finden der richtigen Parameter fur die Entmischung ist dem 

30 Backpropagationalgorithmus ahnlich. 

In dem European Patent application number 93630032.6, date of filing 10. 04. 1992 by the inventors Ehud 
Weinstein and Daniel Yellin wird eine Signaltrennung fiir den Fall einer durch lineare Filter entstandenen 
Mischung unabhangiger und unkorrelierter Quellsignale beschrieben, und es werden zwei Berechnungsmdglich- 
keiten fur die gesuchten Filter gegeben: 1. Eine iterative Losung zweier "standard least square" Algorithmen und 

35 2. ein rekursiver Algorithmus, der im Prinzip dem ersten Equivalent ist, aber weniger Rechenkapazitat bendtigt, 
dafiir aber auch nicht exakt ist Damit wird der von Piatt und Faggin in (J. C. Piatt and F. Faggin: Networks for 
the Separation of Sources that are Superimposed and Delayed, ed. John E. Moody "Advances in Neural Informa- 
tion Processing" Vol.4, Morgan- Kaufman, San Mateo (1992).) vorgeschlagenen Algorithmus begriindet und 
erweitert. In keiner uns bekannten Veroffentlichung wird das Problem einer a priori unbekannten nichtlinearen 

40 Mischung (Mischung durch ein beliebiges nichtlineares System) angeschnitten oder gar gelost 

3. Zugrundeliegendes Problem 

Die Erfindung betrifft ein Verfahren sowie eine Einrichtung zur Signaltrennung eines Signalgemisches der 
45 Form 

+ .- • 

55 M pQty roc 

+ H /„ •••/„ mS U U- ••»«») 

x mx (L-u x )...x mn (l - u n )du x . . . <fu ni 

60 

das aus M unbekannten Quellsignalen xi(t) . . jcm(0 besteht, die mit unbekannten Gewichtsfaktoren g 5 ™,]!,.., j n in 
der angegebenen Weise gemischt sind. Ein Schema dazu stellt Fig. Nr. 1 dar. 1 steht fttr den unbekannten Teil, 
bestehend aus den Quellprozessen, deren Signale durch ein nichtlineares System G 2 gemischt werden. 3 zeigt 
den bekannten Bereich an. In 4 werden aus einer Analyse der MeBdaten die Quellsignale entmischt und stehen 
65 als Ausgange Ui und u 2 zur Verfugung. Das Verfahren kommt fast ohne Annahme uber die Quellsignale und das 
nichtlineare Mischungssystem aus, bis auf die Bedingung, daB eine statistische Analyse der empfangenen Mi- 
schungen aus den Quellsignalen ys(t) genugend Information fur eine Trennung liefert. Beispielsweise wird die 
zeitlich unverzogerte Oberlagerung zweier GauBverteiiter weiBer Rauschsignale m(t) und n^t) mit unterschied- 
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lichen konstanten Faktoren cu und C21 in der Form 

yi(t) = m(t)+ Ci 2 n2(t) 
Y2(t)-n2{t)+ci2ni(t) (1) 

5 

nicht moglich sein, da (yt(t)y2(t + T)) = 0 fur t^O und alle Korrelationen der Form (yVtJy^t + T)) = 0 fur 
m-fn>2. 

Die Beschrankungen werden im Verlauf dieser Beschreibung klarer werden, aber da eine Veranderung der 
g s n,mi t ..,m n durch eine geanderte Aufnahmetechnik moglich ist und da man es in den meisten Fallen mit der 
Trennung von Signalen, die Information enthalten, d h. bei denen eine statistische Analyse hoherer Ordnung 10 
Ergebnisse liefert, wie z. B. eine Trennung von Mischungen aus Sprache und Musik, zu tun hat, werden die 
Beschrankungen meistens ausraumbar bzw. vernachlassigbar sein. 



4. Erfindnung 

Eine beliebige nichtiineare Mischung kann durch die folgende Form beschrieben werden: 

roo r*o 



+ ... 

M 



ZmA* ~ ui)...z mn {t - u n )du x ...du n . (2) 
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wobei die M Quellsignale xi(t) . . . xm(0 unbekannt und mit unbekannten Gewichtsfaktoren gVh,-..,jo in der 
angegebenen Weise gemischt sind und wobei teR ist Dies ist die Form einer Vo 1 terra- Reihe. Entsprechend 
einem Theorem von Frechet kann jedes stetige Funktional mit einer beliebigen Genauigkeit durch eine Volter- 35 
ra- Reihe approximiert werden. Zu dieser Reihe gibt es auch ein zeitdiskontinuierliches Analogon, das eine 
Anwendung der Erf indung auf digitaie Systeme zulaBt 

Im folgenden wird nur der Fall betrachtet, in dem xj(t), ys(t), teR gilt Eine Formulierung mit hdherdimensio- 
nalen Variablen ist moglich, aber der eindimensionale Fall enthalt alle wesentlichen Verf ahrensschritte. 

Das in dieser Erfindung beschriebene Verfahren besteht im wesentlichen aus zwei Schritten: Erstens werden 40 
die durch den wahibaren Grad der Nichtlinearitat bei der Mischung eindeutig bestimmten nichtlinearen Glei- 
chungen ftir ein gleitendes Zeitfenster geldst und die Losungen werden uber die Zeit gemitteit, und zweitens 
wird ein aus einer geniigend groBen Anzahl von unterschiedlichen Kumulanten der geschatzten Ausgangssigna- 
le ui(t) gebildetes Potential minimiert, wobei die zur Berechnung des Potentials notigen Werte aus einem 
gleitenden Zeitfenster einer wahibaren Lange stammeiL 45 

Zu erstens: 

Eine mehrdimensionale Fouriertransformation von 2 ist moglich, und man erhalt: 



Y,(<2) = £gI>,)^( Wi ) + ... 50 

+ H Cyi, .....,„ „ ("•>! 1 • • • » ("1 ) • • • Xm n {•*■ )« 

mi 7ji„ = 1 55 

Genauso ist eine zeitvariable Fouriertransformation mdglich, in der jeweils ein zeitlich begrenzter Ausschnitt 
der GroBe TeR aus dem Gesaratverlauf der Signale herausgeschnitten und fouriertransformiert wird. Dann 
erhalt man: 

60 

Y.&.T) = 53Gij{wi).v i («„r) + ... 

i=l • 

A/ 65 

(W, , . . . , U n )X mi (wi , T) . . . A' m „ (w Ul T). 



3 



BNSOOCID:<DE 19531388C1> 



DE 195 31 388 CI 



Dabei ist vorausgesetzt, daB sich das System G so langsam andert, daB diese Anderung bei der Berechnung der 
gesuchten Mischungsfaktoren vernachlSssigt werden kann. 

Um die Bedingungen zu formulieren, die die Quellsignale erfiillen mussen, damit das Verfahren im ersten 
Schritt angewendet werden kann, werden nun folgende Variablen definiert: 

UJc,No,N ( ,N 2 eN 
coijeR 

(M,...,iNo)e{peN|l<p<M} N ° 
(ji,...jNt)eN N i 
(k l( ..,k N2 )eN N 2 
(a)i ljl ,...<oiNojNi)eR N o N i 

( Ti t j lkl , . . , TiN0jNlkN2)€ R N 0NlN 2 

Die Quellensignale miissen die folgenden Bedingungen erfiillen fur einen genugend groBen Bereich der oben 
definierten Variablen: 

cwn( (3) 
* * 

) =* 0. 

auBer falls 

it — i2 — INo 

und 

<*>iUl ~ • • * — <*>iN0»jNl 

gleichzeitig gilt 

Nun geht es darum, die in der obigen Gleichung stehenden unbekannten GroBen zu berechnen. Man benotigt 
dieselbe Anzahl an Gleichungen wie an unbekannte GroBen. Daher soil mit obiger Gleichung fur ausreichend 
unterschiedlich viele Verzogerungszeiten To bis te (E € N) ein Gleichungssystem aufgestellt werden, d. h.: 

(Ys^~T)Y Sl (^T + T 0 )} = ... 

(Y SQ (^T)Y Si (^T + r E )) = ... (4) 

oder entsprechende Gleichungen, die sich durch das Bilden der Erwartungswerte von beliebigen Produkten der 
zeitvariablen Fouriertransformierten der empfangehen Signale Ys(a>i, . <o n , T) und der zeitverzogerten zeitva- 
riablen Fouriertransformierten der empfangenen Signale Y s (a> l( ..,o>n, T— Tj) ergeben. Dadurch erhalt man ein 
nichtlineares Gleichungssystem, das numerisch gelost werden kann. Naturlich kann man durch die Beriicksichti- 
gung weiterer Werte fQr die Verzogerungszeiten tj mehr Gleichungen als unbedingt notig aufstellen und dann 
die iiberbestimmten Losungen einer Mittelwertbildung unterziehen. 

An einem Beispiel soli nun die Verfahrensweise dieses ersten Schrittes gezeigt werden. 

Das betrachtete System sei derart nichtlinear, daB gilt: 
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GijtwiJ^i.r) 

Gj.i. 3 {cJl. t*l)X\{.txii. T)X2(u.'2, T) 

c;,j,h.L' 2 )Xj(yi.r)x,k. D; 

Oli.a(«i,«=)^i (wi. r)Xa(«4, r) 
G\^ x {uJi.te2)X<i{u:\.T)X\{o;i, T). 



Dann ergeben sich durch das Bilden von Kreuzkorreiationen die folgenden Gleichungen (0<j<E) fiir ©i, 
© 2 €Rund©i^© 2 : 

(Y ffl 2» TJY^cai, ©2, T+ Tj)) = 
ClKwi XXt(® i, T)X,(©,, T + xj)) + 
G ! iXioi XXj(a>i, T)X2(coi, T + Tj)) + 

G'2,1 j(Ol. © 2 )G 2 2,l^(Ci)t, ©2XX,(C0l, T)X|(©t,T + Tj))(X2(©2, T)X2(0)2, T + Tj)) + 

G l 2Ai(a)i, co 2 X> 2 2Ai(<Oi. o>2)(X2(o)i, T)X 2 (©i, T + x^dfsn, T)X,(© 2 , T + Tj)); 

(Y,(©i,©2,T)Yi(©i,© 2 ,T+Tj)) = 

(X l (co,,T)X I (a)i,T+T)) + 

(G l i^©)) 2 ^©,, T)X 2 (<b i ,T + Tj)) + 

(G'2,t,2(©,, © 2 )) 2 (X,(©,, T)X,(6>,, T + Tj))(X2((D2,T + Tj)) + 

(G'2Ai(a>i. ©^RX^©,, T^©,, T + Tj)XX,(© 2 , T)X,(© 2 , T+ Tj)); 

(Y^©!,©*!^©!, ©2, T + Tj)) - 

pC2(<at, T)X2(<»i , T + tj)) + 
(G 2 u (w,)) 2 (X l ((a 1 ,T)X,(a>,, T+ Tj))+ 

(G 2 ^©,, © 2 MX,(©,, T)X,(© 1( T + TjJXX^fflz, TJX^oz. T+ Tj)) + 

(G 2 2Al(©l, ©2)) 2 <X 2 (©1, T)X2(©,, T+ Tj)XX,(©2, T)X,(© 2 , T + Tj)). 

Entsprechend ergeben sich 3 weitere Gleichungen, wenn ©i und © 2 vertauscht werden. 

<Y,(©2,© 1 ,T)Y 2 (©|,©2,T + Tj)) = 
G 2 l,l(©2)(Xi(©2, T)Xi(©2. T + Tj)) + 
G' 1^(©2XX2(©2, T)X2(©2, T + Tj)) + 

G 1 2,l^(©2, ©lJG 2 ^©* ©iXX,(© 2 ,T)X l (©2, T + Tj)XX 2 (© u T)X 2 (©i,T + Tj)) + 

G'lAi^.aiJG 2 ^©* ©iXX 2 (©2,T)X2(©2,T+Tj)XX,(coi,T)X,(©i,T+Tj)); 

(Y,(©2, ©i, T)yi(©2, ©i, T + Tj)) = 
(X,(© 2 ,T)X,(© 2 ,T+Tj)) + 
(G'i^(©2)) 2 (X2(© 2 ,T)X 2 {© 2( T+Tj))+ 

(G'2,lX© 2 . ©t)) 2 ^©* T)X,(©2,T + Tj)XX2(©,, TPC2(©,, T + Tj)) + 

(G'2Ai(© 2 , ©i)) 2 ^©* TJX^Wi T+ Tj)XX,(©i, T)X,(©,, T+ Tj)); 

(Y2(©2, ©l.TJY^©* ©1, T + Tj)) = 
(X2{©2, T)X2{©2, T + Tj)) + 

(G 2 U (© 2 )) 2 (X I (© 2 ,T)X,(© 2 , T+ Tj))+ 

(G 2 2,i^(©2, <»i)RXi(a* T)X,(©2, T+ ^XX^lTJX^©,, T+Tj)) + 
(G 2 2Ai(©2. ©i)) 2 (X2(© 2 , T)X^m, T+ Tj)XX,(©i, T)X,(©,, T + Tj)). 

Fur ein Paar (©t, © 2 ) mit ©i ^ © 2 ergeben sich 6E Gleichungen fOr 1 2 + 4E unbekannte GroBen. 
Ira Fall, da8 ©i, © 2 eR mit ©i =© 2 ist, ergeben sich durch das Bilden von Kreuzkorrelationen die folgenden 
Gleichungen (0 ^ j < E): 

(Y,(©,,©,,T)Y2(©,,©i,T+Tj)) = 
(G 2 i, 1 (©,XX l (©,,T)X I (© 1 ,T+Tj)) + 
(G 1 1^©, XX2(©«, T)X2(©,, T + Tj)) + 

(G>2.i,2(tti, ©iJG 2 ^^©!, ©i)pC t (©i, T)X,(©,. T+^JOX^©,, T)Xz(©,, T+ Tj))+ 
G'2ai(©i» ©iJG^K©!. ©tXX 2 (©i,T)X2(©i, T + Tj)XX,(©,, T)X,(©,, T+Tj)) + 
G'2,1^©,, © t )G 2 2Ai(©!, ©zXX 2 (©i ■ VX4<»u T + tj))(Xi(oi. T)Xi(©i,T + tj)) + 
G'2Ai(©i,©t)G 2 2,i^(©i,©iXX2(©iiT)X2(©i,T+Tj)){Xi(©i,T)Xi(©i,T + Tj)); 
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<Yi(6>i,coi,T)Y,(a>,,G>t,T+Tj)) = 

(X,(a)|,T)X,(a>i,T+Tj)) + 

(G 1 i^(£0,)) 2 <X2(a>,, T)X2(o)i. T+Tj)>+ 

(G'^Uo)!, a),)) 2 {X,(a>,, T)X,(a>i, T+ Tj)XM<»u TJX^o,),, T+ Tj)) + 
(G 1 2 A i(co,.o),)) 2 pC 2 (< 1 >,,T)X 2 (<ai,T+Tj)XX,(o),,T)X,(o>,,T+T j ))+ 
G'^Ucoi, o>,)G'2Ai(a>i, o) l XX|(o>i,T)X,((D,, T+ Tj)XX 2 (o>i, T)X 2 (a>,, T+Tj)>+ 
G'it^ci)!, ©OG'iaiC©!, coiXX2(a),, T)XKw t , T+ Tj)XX,(co,. T)X,((O tl T+ Tj)); 

(Y 2 (o)i, <Oi. TJY^©,, ©i, T+ xj)) = 
(X 2 ({o,.1^X 2 (oj,,T + Tj)) + 
(G*i.i(«)i)XX,(<tt,.-OXi(mi.T+Tj)) + 

(G 2 2,Ua>,, (B,)) 2 ^©,, T)X,(<.),, T+ ^XX^w,, T)X 2 (©, ( T+ Tj)) + 
(G 2 2Ai(«>i, ©i)) 2 <X2(a),, TJX^ca,, T+ Tj)XX,(ca,, T)X,(a>,, T+ Tj)) + 
G 2 ^©,, coOG 2 ^,!^!. a>,XX,((0,, T)X,(co, > T + Tj)XX 2 (<a 1( T)X 2 (a,, T +tj))+ 
G 2 2,i2(a>i,<i)i)G 2 2Ai(<»i, o,XX2(o)i, T^X^O),, T + Tj)XX,(<o,. T)X,(©,, T + Tj)). 

Fur ein Paar (cbi, o 2 ) mit ©i = c*2 ergeben sich 3E Gleichungen fur 6 + 2E unbekannte Gr6Sen. 
Beim Aufstellen dieses Gleichungssystems wurde davon Gebrauch gemacht, daB die Quellsignale die folgen- 
den Bedingungen entsprechend Gleichung (3) erfullen: 



cumiXi (wi , T), X 2 (co7-. T-r)) = 
cum(X l (wt , T),X l ( Wl , T - r), * 2 (a; 2 , T-r)) = 
«m(.Y,(w, , r), .Y,(w x , T - r), A' 2 ("2, T-r)) = 




0 

= 0 
0 

= 0 



fur alle <ou a>2, t e R, und 



cum(Xi(<i)i,T),Xi(a)2»T— x» = 0 
cum(X2(a)t,T) f X2(<02,T— t)) = 0 

fur alle <oi 0)2, x e R mit tot ^ ©2- 
Damit gilt dann: 

(X^tOi^Xstco^TJXtCcouT-tJX^T-T)) = 
(Xi(g>,, T)X!(o)i, T-t)XX2<(02, T)X2(a)2, T- t)>; 
(X2<<Oi,T)X|(a)2,T)X2(a)t,T--t)Xi(oi>2,T-t)) = 
(X2(g>w T)X2(a) I( T-^XXtCtoz, TJX^, T-x)>; 
(X^cot.TpC^cai.TJX^coLT-^X^coi.T-T)) - 
pC2(a)i, TpC^Oi, T- x)XXi(«i f T)Xi(a> lf T-t)>; 

fur alle g>i, o>2€ R. AuBerdem gilt: 

(X^cai.TJX^T^X^uT-^X^coi.T-x)} - 0; 
(X^ca^TJXtCo^TJXiCcOi.T-TjX^i.T-T)) = 0; 

fiir alle a)!, 02 e R mit cot =jfc<02. 

Aufgrund von fehlerhaf ten Aufnahmen (zusatzliches Rauschen) bzw. Fensterungseffekten durch die zeitvari- 
ablen Fouriertransformation mu6 diese Berechnung fiir ein in geeignet groBen Schritten fiber die Signale 
bewegtes Fenster wiederholt werden. Die entsprechenden Ergebnisse werden dann gemittelt, und die gemittel- 
ten Werte nahern sich so den wahren Werten an. Diese Annaherung kann aber viele Daten bendtigen. Eine 
groBe Menge an Daten kann zu Problemen bei einer Realisierung ftihren, daher wird diesem ersten Schritt ein 
zweiter nachgesteilt 

Zuzweitens: 

Dieser zweite Schritt ist eine Anlehnung an die Methoden, die in der Signalverarbeitung mit neuronalen Netzen 
angewendet werden. Man konstruiert sich eine Kostenfunktion, die minimiert wird. Wenn das globale Minimum 
erreicht ist, hat man die optimalen Werte, in diesem Fall der Obertragungsfaktoren, gefunden. Fur die Mimmie- 
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rung einer Kostenfunktion stehen beim heutigen Stand der Technik eine groBe Anzahl von Methoden zur 
Verfugung, z. B. Gradientenabstieg oder Monte-Carlo-Simulationen. Ein Nachteil ailer Verfahren ist aber, daB 
sie in lokaien Minima der Kostenfunktion steckenbleiben konnen; der Vorteil unserer Erfindung ist, daB durch 
die Startwerte aus dem ersten Schritt die Minimierung schon in der Nahe des globalen Minimums beginnen 
kann, denn bei dem Ende der Berechnungen des ersten Schrittes Hegen fur die Obertragungsfaktoren geschatzte 5 
Werte vor, die sich nur wenig von den gesuchten Werten fiir die wahren Obertragungsfaktoren unterscheiden 
und die fiir eine bessere Trennung weiter optimiert werden. 

Mit der Methode der Potentialminimierung werden Effekte von zusatzlichen Rauchprozessen verringert 
Urn die Bedingungen zu formulieren, die die Quellsignale erfullen mussen, darnit das Verfahren im zweiten 
Schritt angewendet werden kann, werden nun folgende Variablen definiert: to 

U.No,Ni€N (5) 
TijeR 

(in . . , i No ) e {p e N| 1 < p < Mfo 
(T il j 1 ,..,r lNO jNi)eR N o N i 

Die Queliensignale mussen die folgenden Bedingungen erfullen fiir einen genugend groBen Bereich der oben 
definierten Variablen: 

20 

cum( (6) 
Xjj (i — 77, j x ),..., x^ (t — ), 

• rt *.v 0 t *n 0 Ji )> • * • 1 ~ r tV 0 ,jVj ) 

) = o } 

auBer falls 
ii = i 2 = . . . = i NQ 

gilt 35 

Fur den hier betrachteten Fall, die richtigen Obertragungsfaktoren zu finden, bilden wir eine Kostenfunktion 
derart, daB die Ausgange eines adaptiven Filters ui(t) dieselben Bedingungen erfullen mussen wie die Quellsigna- 
le *i(t), d. h. fOr eine gentigend groBe Anzahl von Werten der in 5 definierten Variablen werden die Quadrate der 
Kumulanten aus 6, dabei werden die Quellsignale x durch die Ausgange des Filters u ersetzt, aufsummiert, so daB 
eine Kostenfunktion entsteht, die minimiert wird Optimiert werden dabei die Werte der Obertragungsfaktoren. 40 
Nach jedem Optimierungsschritt mussen die Werte fiir die Ui(t) neu berechnet werden. Diese Berechnung kann 
im Zeit- oder Frequenzbereich erfolgen. 
Zum Beispie! werden fur den Fall, in dem folgende Mischung vorliegt: 

Y t (a>) - X^oJ + Gi^X^coJ-M^coJX^caJX^Q)); 45 
Y 2 (a>) = X 2 (ci)) + G2 I (o>)Xi(a)) + K(a>)Xi(e))X2(c t )) 

die Ausgange mit den geschatzten Werten fiir die Obertragungsfaktoren G*i2» G* 2! und K* berechnet durch 
Auflosen des Gleichungssystems: 



25 



30 



Y^a)) = U^+G^iC^JU^coJ + K^^UiCcoJU^o)) 



50 



nach Ui(a>) und U^a)). 

Die durch das adaptive Filter zu minimierende Kostenfunktion V kann mit den Ausgangen ui(t) und U2{t) 55 
gebildet werden: 

V= [cutn{u l (i)^(t + T))] 2 + 

r<zB 60 
[cUTT7.(u l (t),U l (t + T),U 2 (t + r))] 2 + 

(cum(w 2 (f), u 2 {t + r), Ul (t + r))] 2 + 

[cum(u x (t) t u 2 (t), u^t + r), u 2 (t + r))] 2 + 65 
[cur>i(u 2 (t),u x (t),u 2 (t + r)\uj(i + r))] 2 , 
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wobei B eine endliche Menge reeller Zahlen ist 

5. Gewerbliche Anwendbarkeit 



Die durchzufiihrenden Berechnungen bendtigen eine hohe Rechenkapazitat, die aber schon heute mit Mikro- 
prozessoren bereitsteht. Eine in parallelen Schritten ausgefuhrte Verwendung des angegebene Verfahrens 
erdffnet auch fur hdherfrequente Probleme ein Echtzeitrealisierung. Im offline Betrieb laBt sich das Verfahren 
auf frei programmierbaren Mikroprozessoren impiementieren. 

Anwendungsmdglichkeiten liegen in einer Signaltrennung als Vorverarbeitungsschritt fur eine automatische 
Spracherkennung, bei der Konstruktion technischer Hdrhilfen, sowie in alien anderen akustischen Messungen, 
einschlieBlich die der Unterwasserakustik. 

AuBerdem sind medizinische Anwendungsmdglichkeiten bei Mehrkanalmessungen wie EEG oder MEG 
(Fig. Nr. 2) moglich, sowie bei Antennenfeldern fur die Aufnahme elektromagnetischer Strahlung. Wie in der 
Fig. Nr. 2 dargestellt werden magnetische Dipole im Gehirn 2 durch SQUIDs 3 gemessen, und die Signale 
werden einer Auswerteeinheit 1 zugefiihrt 

Diesen und alien weiteren denkbaren Anwendungsmdglichkeiten ist gemeinsam, daB sie der bisher linearen 
Theorie verschlossen bleiben, weil eine lineare Signaltrennung zu unbefriedigenden Ergebnissen fuhrt, denn 
aufgrund der Detektortechniken haben fast aile Empfanger eine nichtlineare Kennlinie. Weiterhin sind bei einer 
Formulierung mit mehrdimensionalen Variablen Anwendungen in der Optik bzw. Bildverarbeitung moglich, in 
der auch bei der Aufnahme nichtlineare Effekte auftreten. 

Daruber hinaus konnen die berechneten Werte auch zu einer Analyse des mischenden Systems benutzt 
werden, z. B. zur Quellenlokalisierung oder Bewegungsdetektion. 

6. Vorteilhafte Auswirkungen 

Neben den schon in dem Abschnitt Gewerbliche Anwendbarkeit genannten mdglichen Erneuerungen gibt es 
noch eine weitere Verbesserung: Urn mit der bisher linearen Methode zu arbeiten, werden Empfindlichkeitsbe- 
reiche der Detektoren auf ihren linearen Bereich verkieinert Dies ist bei der nichtlinearen Methode nicht notig. 



7. Ausfuhrung 

Stellvertretend fur alle anderen Probleme wird der Fall der Trennung akustischer Signale behandelt Entspre- 
chend Fig. Nr. 1 stellt der mit 1 bezeichnete Block dann einen unbekannten ProzeB bestehend aus zwei Schall- 
quellen dar, deren Signale durch den Raum zu zwei Mikrofonen gelangen, deren aufgenommene Signale 
aufgezeichnet werden und die Grundlage der der Erfindung entsprechenden Berechnungen bilden. Der Block, 
der mit 3 gekennzeichnet ist, enthalt alle bekannten GrdBen. Das System G ist nichtlinear, weil zum einen die 
Raumiibertragung nichtlinear sein kann, zum anderen die Mikrofone mit ihren Verstarkern eine nichtlineare 
Kennlinie besitzen. Die berechneten Signale ui und u 2 sind den Quellsignalen k x und x 2 ahnlich, so daB ein 
automatisches Spracherkennungssystem mit der aus der Mischung herausgetrennten Sprache funktioniert Dies 
ist in der Fig. Nr. 3 sichtbar. Dort sind in den oberen beiden Kasten die Mischsignale von Sprache und Musik 
aufgezeigt wie sie sich bei einer Mikrofonaufnahme in einem kleinen Raum ergeben. Im untersten Kasten ist das 
aus dem Gemisch herausgetrennte Sprachsignal am Ausgang der Trennungseinrichtung, dessen Berechnung auf 
einem der Erfindung entsprechenden Verfahren basiert, dargestellt; mit dem untersten Signal ist eine fehlerlose 
automatische Erkennung mit einem kommerziellen automatischen Spracherkenner moglich, wahrend derselbe 
Erkenner bei den oberen beiden Signalen jeweils drei Fehler liefert 



8. Figurenbeschreibung 

Fig. Nr. 1 : 1 steht fur den unbekannten Bereich, bestehend aus den Quellprozessen, deren Signale durch ein 
nichtlineares System G 2 gemischt werden. 3 zeigt den bekannten Bereich an. In 4 werden aus einer Analyse der 
MeBdaten die Quellsignale entmischt und stehen als Ausgange ui und u 2 zur Verfugung. 

Fig. Nr. 2: Magnetische Dipole im Gehirn 2 werden durch SQUIDs 3 gemessene, und die Signale werden einer 

Auswerteeinheit 1 zugefiihrt . . . . 

Fig. Nr. 3: In den oberen beiden Kasten sind die Mischsignale von Sprache und Musik aufgezeigt, wie sie sich 
bei einer Mikrofonaufnahme in einem kleinen Raum ergeben. Im untersten Kasten ist das aus dem Gemisch 
herausgetrennte Sprachsignal am Ausgang der Trennungseinrichtung, dessen Berechnung auf einem der Erfin- 
dung entsprechenden Verfahren basiert, dargestellt; mit dem untersten Signal ist eine fehlerlose automatische 
Erkennung mit einem kommerziellen automatischen Spracherkenner moglich, wahrend derselbe Erkenner bei 
den oberen beiden Signalen jeweils drei Fehler liefert 

Patentanspriiche 
1. Verfahren zur Separierung eines Signalgemisches der Form 
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+ ... 



roo roo 

+ Z, / **7n ^ m«(«I> •■•.«») 
Zm t (t - Ui) . . - X mn (i - U n )<^l . . . du n , 



Hierzu 2 Seite(n) Zeichnungen 
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das aus M unbekannten Quellsignalen xi(t) . . . xmM besteht, die mit unbekannten Gewichtsfaktoren is 
S\lu • in * n der angegebenen Weise gemischt sind, fixr die Variable t soli gelten: teR, dadurch gekenn- 
zeichnet, daB ersten die durch den Grad der Nichtlinearitat bei der Mischung eindeutig bestimmten 
nichtlinearen Gleichungen, die durch das Bilden von einer geniigend groBen Anzahl von unterschiedlichen 
Erwartungswerten von Produkten der zeitvariablen Fouriertransformierten der empfangenen Signale 
Y s (a>i,.., a>n, T) und der urn Tj zeitverzogerten zeitvariablen Fouriertransformierten der empfangenen 20 
Signale Y s (wi, . . ., co 0 » T—tj) entstehen, fur ein gleitendes Zeitfenster gelost werden und die Losungen tiber 
die Zeit gemittelt werden und daB zweitens ein aus einer geniigend groBen Anzahl von unterschiedlichen 
Kumulanten aus den geschatzten Ausgangssignalen ui(t) und aus den um Tj zeitverzogerten geschatzten 
Ausgangssignalen ui(t— tj) gebildetes Potential minimiert wird, wobei die zur Berechnung des Potentials 
notigen Werte aus einem gleitenden Zeitfenster einer wahlbaren GroBe stammen. 25 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB uber die Werte der Nichtlinearitat g irgendwel- 
che a-priori-Informationen bestehen, so daB die Berechnungen vereinfacht werden konnen. 

3. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB M = 2 ist. 

4. Verfahren nach Anspruch I bis 3, dadurch gekennzeichnet, daB uber die Werte der Quellsignale % 
irgendwelche a-priori-Informationen bestehen, so daB die Berechnungen vereinfacht werden konnen. 30 

5. Einrichtungen, dadurch gekennzeichnet, daB sie die Verfahren der Anspriiche 1 bis 4 in digitaler oder 
analoger Weise ftir eine Signaltrennung benutzen. 
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